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SUMMARY 


Owing to the quantity of material to be presented, this 
report is divided into several volumes. Volume 1 introduced 
the series, described the analytical procedures in detail, and 
presented and discussed the results for Shot 10. Volumes 2 
and 3 presented and discussed the results for Shots 9 and 8. 
This volume presents and discusses the results for Shot ll. 

A subsequent volume will compare the results of the four 
experiments. The method of analysis is common to all four 
experiments and is described in detail in Volume 1 only. 

So that the results from the four experiments may be 
easily compared, they have been scaled to remove the effects 
of varying atmospheric conditions. (Results are scaled to a 
1 kg charge weight and a standard atmosphere of dry air at 
15°C at sea level.) For the most part, only scaled results 
are presented. Exceptions include some derived pressure-time 
histories, which are compared to actual gauge measurements 
made in the experiment. 

Results are presented in SI units, even though the 
experiments were originally laid out in British units. Only 
distance and time measurements are affected, however, as 


velocity, density, and pressure results are presented as 


dimensionless ratios. A distance units conversion scale is 
: : included to convert between SI units (meters scaled to a 


: 1 kg charge) and British units (feet scaled to a 1 1b charge), 


plus a time scale factor and scale factors to convert pressure 
ratios to both British and SI pressure units. Scale factors 
which may be used to compute the distance and time values 


actually observed under the ambient conditions of each shot 


are also provided. Dimensional pressure units are used for 


the results presented at gauge locations. 
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CHAPTER 1, SHOT 11 ANALYSIS 


1.1 Introduction 


This is the fourth volume in a series which presents the 
particle trajectory analysis results from four experiments 
(Dipole West Shots 8, 9, 10 and 11) carried out to obtain 
information on the interaction of spherical blast waves with 
real and ideal reflecting surfaces. A general description of 


the project can be found in Volume 1. The results presented 


in this volume are for Shot ll. 


In each experiment, photogrammetrical studies were made 


of the shock fronts (refractive image analysis, RIA), and of 
the motions of smoke puff particle tracers (particle trajectory 
analysis, PTA). The refractive image analysis results were 
reported by Dewey et al. (1975) and results of the particle 
trajectory analysis are presented in this report. The method 
of particle trajectory analysis, common to all four shots is 


described in detail in Volume 1 only. 


1.2 Description of Shot Il 


Dipole West Shot 1llwas fired on November 8th, 1973 by 


the Ballistics Research Laboratories at the Defence Research 


Establishment Suffield, in Alberta, Canada. Two 1080 lb (490 kg) 
spheres of Pentolite were detonated simultaneously, to within 


5 microseconds, at nominal charge heights of 25 and 75ft 


(7.6 and 22.9m) over a relatively rough ground surface. 


Particle trajectory data were gathered by photographing 
the movement of smoke puffs formed in a vertical plane running 
out from ground zero at 6.7° south of west. A WF5 camera 
Operating at about 3500 frames per second was positioned 
30ft (9.2m) above ground level at a position 600ft (183m) 
due south of ground zero (GZ), the point on the ground 
vertically beneath the charges. 

Table 1 gives the field survey data for the event, and 
Figure 1 shows a plan view of the layout. The dashed line 
represents the approximate line of sight of the WF5 camera. 
Figure 2 shows the field of view of this camera. 

The smoke puff grid was made up of 20 columns of 12 puffs 
each, hung vertically on strings. The vertical spacing of 
puffs was 5ft, beginning 3ft above ground level and ending 
at a height of 58ft. The horizontal spacing of the columns 
of puffs was 10, 7 or 5ft, depending on the distance from 
ground zero, beginning at about 25ft and ending at about 140ft 
from GZ. Of the possible 240 smoke puffs, 232 detonated 
successfully. A good film record was obtained, except that 
several of the white smoke puffs in the bottom row could not 
be very easily distinguished against the background of snow- 
covered, ploughed earth. 

This report describes the analysis of the smoke puff 
data collected for Shot 11, and presents and discusses some 


of the results of that analysis. 
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1.3 Camera calibration and data reduction 


The calculated camera position coordinates and orientation 
angles for Shot 11 are presented in Table 2, together with 
the positions of photomarkers transformed from one frame of 
the film just before detonation to an object plane defined 
aS passing through ground zero and being normal to the camera 
orientation axis. The differences ("Shifts") between the 
object plane positions of the transformed calibration points 
and their positions computed from the field survey data are 
given in Table 2. The object plane positions of the calibra- 
tion points computed in these two ways are also shown in 
Figure 3. 

The camera calibration procedure, described in detail in 
Volume 1, ensured that selected photomarker images (Pl to 
P5) transformed to the object plane in a way which matched 
them exactly to the positions computed using the survey data. 
These reference photomarkers for Shot 11 are indicated in 
Figure 3 using large circles: namely, Pl = Wl, P2 = W3, 
and P3 = 300Wl. The separation distance between P4 = P3 = 
300Wl and P5 = 300W2 was also used as a calibration parameter. 
The probable reason for the shifts seen for photomarkers VP1A 
and VP1B was discussed in Volume l. 


The image positions of two reference photomarkers 


(VP3B and 300W2) and all smoke puffs were measured frame-by- 


frame over a time interval corresponding to the approximate 


ll 


duration of the positive phase of the blast waves (film frames 
10 to 375), and were transformed to distances in the object 
plane by matching the reference marker positions to their 
positions transformed from the calibration frame. These data 
were again transformed from the object plane to the smoke puff 
plane which was assumed to pass through "corrected" ground zero; 
to be vertical, and to run 6.7° south of west from GZ. 

The x-y coordinate system in the smoke puff plane was 
the same for Shot 1l as for the other shots, except that the 
corrected value for ground zero was displaced 0.7ft from the 
surveyed ground zero, in a direction approximately 34° south 
of west. The corrected ground zero was defined to have the 
same elevation as the surveyed ground zero, but was located 
directly under the midway point between the two charge centers. 
As for the other shots, all data in the output plane are 
plotted with the x coordinate reflected, i.e. with positive 
values of x to the right hand side, as if the smoke grid had 
run to the right of the charges rather than to the left as 
seen in the film images. 

A time was assigned to each film frame using the 1 ms 
timing marks placed on the film during its exposure. The 
film timing method was described in Volume 1, and the 
complete set of film timing data used for Shot 11 is provided 
in Table 3. 

Figure 4 shows the positions of the detonated smoke puffs 


at a time prior to the detonation of the two charges. These 
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positions are in the plane of the charges and the smoke puff 
grid, as described above. The smoke puff plane was not 
exactly parallel to the camera image and object planes 
(Figures 2and 3), and various geometrical corrections were 
applied to make the transformation between them. The puffs 
enclosed in parenthesis were not visible in the earlier film 
frames because they were concealed by the photomarkers, 


but were seen later. The puffs which are underlined in the 


figure were seen in their initial positions only - no 

trajectory data could be obtained. Charge positions in the 
figures are plotted as if they were positioned exactly above 
the corrected ground zero origin. The data shown in Figure 


4 have not been scaled. 


1.4 Data scaling and trajectory fitting 


The position-time histories of individual smoke puffs 
were extracted from the frame-by-frame positions of the smoke 
puff grid, and then scaled to standard atmospheric conditions 
and charge weight. A change to SI units was made at this 
point in the analysis. The resulting trajectories were edited, 
and then smoothed by fitting polynomial functions. 

Particle trajectory data were scaled by dividing all 
distances by Sachs scaling factor S = Y(WPo)7 (WP) and 
multiplying all times by the factor C/(C5S), where C is the 


ambient sound speed computed for Shot 11. Data used to compute 
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C and S, and define the scaled event, are listed 


the computed values of C and S. 


Ambient temperature, T 
Ambient pressure, P 
Relative humidity, RH 
Computed vapour pressure, VP 
Ambient sound speed (eo 
Charge weight, W 
Sachs scaling factor fs) 
Standard charge weight, WwW 
Standard pressure, P 
Standard temperature, T 
Standard sound speed, Cc 
(dry air) 


-19.11 °c 
94.34 kPa 
60.0% 

0.08 kPa 
319.639 m/s 
489.9 kg 


8.0730 


1.0 kg 
101.325 kPa 
US. 26 


340.292 m/s 


below with 


(-2.4 °F) 


(13.683 PSI) 


(0.6 mm Hg) 
(1049 ft/s) 
(1080 lbs) 


(2.2 lbs) 
(14.7 PSI) 
(59 °F) 


(1116 ft/s) 


The results presented in this report therefore apply to 


a scaled event which is the detonation of two 1 kg charges in 


a standard atmosphere. The scaled heights of burst for Shot 11 


were 0.905 m and 2.797 m, and the scaled charge separation 


divided by two, was 0.946 m. 


Figure 5 shows the scaled particle trajectory data for 


Shot 11 in the smoke puff plane with positions measured 


horizontally and vertically from 
mately 27,247 puff positions are 
the raw trajectory data have not 


The raw particle trajectory 


obvious data processing errors, 
14 


corrected ground zero. 
represented. 


been smoothed. 


As represented, 


data were edited to remove 


such as a single point widely 


Approxi- 


ih 


displaced from its trajectory for one or two frames. The 
trajectory of each puff in turn was then smoothed by least 
squares fitting simple polynomial expressions separately to 

both the x and y coordinate data, these being discrete 

functions of frame time. The adequacy of each fit was determined 
by examining on the same graphical output, plots of both the 

raw trajectory data and the fitted curve. For Shot 11 this 

Meant examining and adjusting 464 such plots, at least two 

or three times each. 

For a given puff, the first step in fitting the raw 
trajectory data was to set the time of arrival of the shock 
front first hitting the puff. The data at subsequent times 
were fitted with polynomial functions, as described in Volume 
1, paragraph 2.5. The first derivatives of the fitted 
functions were also calculated at a series of times for use 


in later calculations of particle velocity. 


1.5 Regionalization and shock strength calculations 


Five regions were defined in the smoke puff plane on 
the basis of the shock front which first struck the puffs in 
a particular region. These are shown in Figure 6. The 
regions were bounded by the triple point trajectories 
measured using refractive image analysis (Dewey et al., 
1975). Regions 1 and 2 are those in which the smoke puffs 
were first hit by a spherical primary shock front, and regions 


3, 4, and 5 are those in which the puffs were first hit by a 


1S 


Mach stem. 
In each of the five regions, the shock trajectory data 
obtained from the first movement of the smoke puffs were 
fitted to a function of the form 


j r(t) =A + Bt + C log (1 + t) + DVlog (1 + t€) , 


where r is the shock radius, t is the time after detonation, 


and A, B, C and D are the fitted coefficients. The shock 


velocities were calculated by differentiating this function. 


The peak particle velocity, Ve , peak density, De , and peak 


{ hydrostatic overpressure, Ps , aS functions of shock radius in 


each of the five regions, were calculated from the shock 


velocity using extensions of the Rankine-Hugoniot equation. 
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Details of the shock radius calculations etc. are described 


in Volume 1, paragraph 2.6. 


1.6 Particle velocity calculations 


Particle velocities were computed using the methods 


described in Volume 1, paragraph 2.7. 


1.7 Density and hydrostatic overpressure calculations 


Densities and hydrostatic overpressures in the smoke puff 
plane were calculated by the method described in Volume l, 


paragraph 2.8. Results in both cases represent average values 


over cells defined by four adjacent smoke puffs. 


1.8 Surface representation 


Surfaces were fitted to the times of shock front arrival 
and to the fields of particle velocity, density and hydrostatic 
overpressure at a sequence of times. All data were interpolated 
onto a common regular Euleurian grid. Fields of dynamic 
pressure were computed from surface-interpolated particle 
velocity and density results. Contour plots were generated 
for all surfaces at selected times, and time histories computed 
at several fixed locations. The methods used were identical 


to those described for Shot 10 in Volume 1, Chapter 3. 


1.9 Pressure and total pressure time-histories 


To permit a direct comparison between results obtained 
from the particle trajectory analysis and measurements made 


using side-on and face-on pressure transducers, the hydro- 


static and total overpressure time-histories were calculated 
at those locations coincident with gauge positions within the 
smoke puff grid. Dynamic pressures and hydrostatic over- 
pressures obtained from the particle trajectory analysis were 
used to compute the total pressures after application of a 
compressibility correction. This correction is a function of 
the local Mach number and its form depends on whether the Mach 
number was greater or less than unity. The time varying 
hydrostatic and total overpressure impulses, determined by 
integrating the pressure time histories, were also calculated 
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and compared with similar integrations of the electronic 


transducer data. 


The methods used to calculate the total pressures and 


See wees, eRe 


the impulses are described in detail in the addendum to 
Volumes 1 and 2, which is incorporated in Volume 3. In 
cases where the leading edge of a time-history curve was 


rounded, impulse integrations were done using data inter- 
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polated linearly between the peak parameter value determined 
at the time of arrival, and a point on the time-history curve 
subsequent to the time of arrival. The second of these two 
points was chosen in a manner which ensured a minimum 
difference in slope between the interpolated and computed : 


sections of the time-history data. 
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CHAPTER 2. SHOT 11 RESULTS 


2.1 Times of shock front arrival 


The measured initial puff positions, the times of first 
shock arrival, and the peak particle velocities obtained by 
differentiating the functions fitted to the particle 
trajectories are presented in Table 4. Puff position is 
given relative to corrected ground zero as origin, with hori- 
zontal and vertical axes. Puff position and the time of 
arrival of the first shock are given both as observed and 
scaled. Particle velocities listed are derivatives of the 
fitted puff trajectories at the times of shock arrival, and 
are expressed in Mach units. Expressed this way, the particle 
velocities are the same scaled as unscaled. Also listed are 
the initial radial puff positions (scaled) and region codes. 

Shock front data determined from the first movement of 


the smoke puffs, i.e. calculated from the time-of-arrival 


data in Table 4, are listed in Tables 5.1 - 5.5. Each table 
corresponds to one of the 5 regions used. Listed are the 
observed and fitted unscaled shock trajectory data, the 
scaled fitted shock trajectory data, and the computed shock 
velocities and peak parameters associated with shock strength: 


peak hydrostatic overpressure in atmospheres and kilopascals, peak 
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particle velocities in Mach units, and peak density ratios. 
Given as ratios, these peak parameters are the same scaled 
as unscaled. Pressure given in kilopascals in the tables 
refers to the unscaled (observed) case only. 

The shock front radius versus time data derived using 
particle trajectory analysis (PTA) are also shown in Figures 
7.1-7.3 for the two primary fronts, the two Mach stems at 
the interaction plane, and the ground Mach stem, respectively. 
They are compared to corresponding data derived from refractive 
image analysis (RIA) reported by Dewey et al. (1975). The 
refractive image analysis results were obtained using photo- 
grammetry against a striped canvas backdrop and they describe 
the shock as it travelled in a direction almost diametrically 


opposite to the direction of the smoke puff grid. 


2.2 Shock strengths 


Peak particle velocities calculated from shock front 
velocities are shown in Figures 8.1- 8.3 for the primary 
fronts, interaction Mach stems, and the ground Mach stem. 

This method of determining peak particle velocities was 
labelled method 1, and the data plotted correspond to 

those listed in Tables 5.1-5.5. The results in the figures 
are compared with those previously obtained using refractive 
image analysis (RIA). In the case of the primary shock fronts, 


results are also compared to those of Brode (1957) for TNT. 
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In Volume 1 other methods of determining shock strengths 
in the various regions were described. It was demonstrated 
that method 1 was clearly the most accurate, and in the present 
volume shock strengths calculated using methods 2 and 3 are 
not reported. For this reason Figures 9, 10 and 11 and 


Table 6 do not appear in this volume. 
2.3 Particle velocity fields 


The calculated particle velocities in the plane of the 
smoke puffs are shown as vectors in Figures 12.1 through 12.8, 
for various times after the detonation. All times and 
positions are scaled to a 1 kg charge in a standard atmosphere. 
The particle velocity vectors represent the derivatives of 
the smoothed particle trajectories, and their magnitudes may 
be judged using the standard vector shown on each figure. All 
velocities are measured in Mach units, relative to the standard 
sound speed. Puffs not yet struck by a shock wave are 
represented by small circles (zero velocity). 

Numerical data corresponding to Figures 12.1 - 12.8 are 
listed in Tables 7.1 through 7.11, along with scaled radial 
positions of the puffs, and region codes as defined in Figure 6. 
Conversion factors are given at the foot of each table, which 
may be used to convert the scaled data in the tables and figures 


back to their original unscaled values. 
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2.4 Density and hydrostatic overpressure fields 


Calculated average relative densities throughout the 
smoke puff plane are depicted graphically in Figures 13.1 - 13.4, 
for various times after the detonation. All time values are 
scaled. Cell positions are scaled and are given relative to 
the corrected ground zero as origin with horizontal and 
vertical axes. The calculated densities may be judged using 
the density shading scale shown on each figure. Density is 
given as a ratio, relative to ambient density. Cells not yet 
struck by a shock wave and cells in which the density has 
dropped to a value less than ambient density are shown blank. 

Corresponding numerical data are listed in Tables 8.1 - 
8.8 along with radial cell positions computed according to 
the regions defined previously. Numerical data describing the 
fields of hydrostatic overpressure are similarly listed in 
Tables 9.1 - 9.8. The pressure results for a given cell were 
obtained by multiplying the density results for that cell by 
a factor determined by the strength of the shock which first 
traversed the cell and which then remained constant, i.e. by 


assuming isentropic flow after the first shock. 
263 Times~of-arrival surface 


Figure 14 shows a perspective view of the surface fitted 
to the original smoke puff positions and the observed times of 


first shock front arrival, i.e., to the data listed in Table 4. 
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The grid mesh size is 0.1 by 0.1 meters (scaled), about 

2.5 feet square (unscaled), or about half that of the original 
smoke puff grid. The charge positions are indicated on the 
vertical distance axis. 

The times-of-arrival surface is smooth enough to permit 
contouring, the contours in this case (isochrones) representing 
shock front shapes at different times, as shown in Figure 15, 
but the surface is not smooth enough to permit the calculation 
of gradient vectors which could be used to compute shock 
velocity vectors and shock strengths over the new grid. 

Two attempts were made to obtain contours of shock 
strength. In the first, the times-of-arrival surface was 
smoothed by least-squares fitting low-order, one-dimensional. 
polynomial functions to the time-of-arrival data along each 
grid row and column separately, and computing the derivatives 
of the fitted functions to obtain the associated components 
of the surface gradient vectors. Shock velocity vectors were 
obtained from the time-of-arrival gradients, and from these 
peak particle velocities were computed. The peak particle 
velocity (shock strength) surface is shown in Figure 16. 

The contours of this surface (not shown) did not exhibit any 


discontinuities across the boundaries of the shock front 


regions, as they would if surfaces had been fitted to the time 
of arrival in each region separately. 
The results of a second method used to compute shock 


strength contours are shown in Figure 17. These were obtained 


eo 
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by interpolating shock radius at each value of peak particle 


velocity shown, for each shock front region in turn, using 

the peak particle velocity versus radius curves shown in 
Figures 8.1-8.3. Arcs of circles with these radii, centered 
on the appropriate points along the vertical charge axis, were 
then drawn in the regions to represent shock strength contours. 
These peak value contours are discontinuous across triple 
point locii and other region boundaries. As a result, some 
horizontal lines are crossed twice by the same contour or, in 
other words, indentical shock strengths can be found at two 


locations the same vertical distance from a reflecting surface, 


but at different radial distances from the vertical charge axis. 


2.6 Field surface contours 


Contours of equal particle velocity, density, hydro- 
static overpressure, and dynamic pressure in the blast waves 
were determined for a series of times, using surfaces fitted 
to the various measured data fields at those times. Sample 
results are shown in Figures 18 through 21 at scaled times 
of 2.5, 4.0 and 9.0 ms. The shock fronts shown in these 
figures are obtained from the time-of-arrival surface (as were 
those in Figure 15). Field contours such as those shown can 
be drawn for any scaled time between 0.5 ms and 12.7 ms. 


It should be re-stated that all of these results were 


obtained from the photography of the smoke puffs only and do 
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not rely on the results obtained using the refractive image 


analysis (Dewey et al., 1975). 
2.7 Time histories 


By mapping the physical properties of the blast waves 
at short time intervals its was possible to determine the time 
histories of these properites at any selected fixed position 
within the smoke puff grid. This was done at 15 fixed 
locations, three in the primary region of the lower charge and 
four in each of the three Mach stem regions, as shown in 
Figure 22. At each distance from the vertical axis through 
the charges in the Mach stem regions, each of the time 
history stations is approximately the same distance from either 
the interaction plane or the ground plane. (Particle velocity 
time histories could be interpolated closest to the grid edges 
because these were measured at puff locations, whereas the 
density and pressure data were measured at cell centers). 

Time histories of particle velocity, density, hydrostatic 
overpressure and dynamic pressure at these locations are given 
in Figures 23 to 26. Time-histories of these physical properties 
of the blast wave can be provided at any other location within 
the smoke puff grid, on request. 

The vertical line which forms the leading edge of a time- 
history plot represents the interpolated time-of-arrival of 
the first shock at the given location, and the height of this line 


represents the peak parameter value derived from the shock velocity 


analysis. 25 


The dynamic pressures plotted in Figure 26 are maximum 
values, computed using both the x and y component of particle 
velocity. Similar plots were made of the horizontal 
components of dynamic pressure, but the differences were not 
significant since the y components of particle velocity at 
these locations were small. Other locations could have been 
chosen at which the y components would not have been insigni- 
ficant. 


Time histories for hydrostatic overpressure and total 


pressure are also plotted in Figures 27.1 to 27.4 for stations 


at the nominal positions of field-mounted pressure gauges on 
the "60 foot gun barrel". The gauges on this gun barrel were 
mounted at nominal elevations of 10, 20, 30, 40, 47, 50 and 
53 feet. The time histories at these locations are compared 
to the gauge measurements (Keefer and Reisler, 1975). The 
total pressures were calculated in the manner described in 
the addendum to Volumes 1 and 2 which is incorporated in 
Volume 3. The variation with time of the integrated pressure 
(pressure impulse) is also shown in these figures, compared 
with similar integrals of the gauge data (Keefer and Reisler, 


1975). 
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CHAPTER 3, DISCUSSION 


3.1 Particle trajectory analysis, Shot 11 


The methods used to analyze the smoke puff trajectories 
on Shot 11 were identical to those used for Shots 8, 9 and 
10 and described in detail in the first three volumes of 
this report. The results for Shot 10 clearly indicated the 
superiority of one of several methods of analyzing shock 
strength, and only the results of this method have been 


reported for Shots 8, 9 and ll. 


3.2 Primary shock strength of upper charge 


The refractive image analysis of the shock fronts 
described by Dewey et al (1975) did not provide any information 
about the primary spherical shocks from the upper charges, and 
it was assumed that these charges had detonated satisfactorily. 
This assumption was validated for Shots 10 and 9 by the 
analysis of the particle trajectory time-of-arrival measurements. 
In Figure 7.1 the primary shock radii are plotted versus time 
for the upper and lower charges of Shot 11, together with 
the results obtained for the lower charge by the refractive 
image analysis. All three curves appear to be identical. 

The limited range of the data obtained for the primary shock 
from the upper charge did not permit an accurate calculation 


of the variation of the shock strength with distance. 


However, the results for the lower charge, in Figure 8.l, 
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show a similar shock strength variation with distance to that 
obtained from the refractive image analysis and one which is 


very similar to Brode's (1957) calculation for TNT. 


3.3 Comparison of Mach shock strength over different surfaces ] 


The refractive image analysis of Shot 11 (Dewey et al., 
1975) showed what appeared to be a significant difference 
between the strengths of the Mach shocks over the rough 
ground and beneath the interaction plane between the two 
charges. The results of the particle trajectory analysis 
given in Figures 8.2 and 8.3 show a similar but smaller 
difference. The RIA measurements were made as close as 
possible to the reflecting surfaces, 0.5 m above the ground 
plane and 0.2 m below the interaction plane, whereas in the 
PTA case the results represent an average of measurements made 
at puff positions at heights ranging between 1.0 and 7.0 m. 
Obtaining particle trajectory data near the ground surface 
was a particular problem on Shot 11 since 13 of the 20 smoke 
puffs in the bottom row either failed to detonate or could not 
be distinguished in the film images (they were white smoke 
puffs placed against a background of the snow-covered, 
ploughed earth). The results in Figures 8.2 and 8.3 therefore 
indicate that the difference in shock strength over the 
ground compared with that at the interaction plane may be 
dependent on the height above the ground - not an unexpected 


result. 
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In addition, determination of Mach shock strength from 
measurements of the times of shock arrival at smoke puffs at 
various distances from the reflecting surfaces is difficult 
because an assumption must be made about the exact shape of 
the Mach shock front, in order to correctly assign shock 
radius values at smoke puff positions. At or near a reflecting 
surface the problem of shock shape is not so important as it 
is assumed that the shock is perpendicular to the surface. 
Details of this aspect in the PTA case and the manner in which 
the problem was dealt with for Shots 8 through 1l are described 
in Volume 1 of this report. 

The possible magnitude of the difference between Mach 
stem shock strength over the rough ground and at the interaction 
plane is best illustrated in Figure 17, which shows peak 
particle velocity isotachs derived from the RIA results 
presented in Figures 8.1- 8.3. It can be seen, for example, 
that the Mach 0.5 isotach at the interaction plane is signi- 
ficantly farther from the vertical charge axis than it is at 
the ground plane. The apparent discontinuity in the Mach 
0.5 isotach at the plane y = 0.8 cannot be maintained and 
there will be a transfer of energy between the two shock 
waves causing the upper Mach stem to slow down and the lower 


to strengthen. 
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3.4 Resolution of time histories 


The time histories of density and pressure given 
in Figures 24 and 25 do not always show a sharp rise at 
the shock front. This is not a real effect but one inherent 
to the present method of particle trajectory analysis, which 
does not permit a high resolution of density in space or in 
time because the average density of the air within a rectangular 
cell defined by four smoke puffs cannot be calculated accurately 
until the shock has completely traversed the cell. The 
time of complete traversal may be as much as 5 ms, or 0.6 ms 
when scaled. 

For the same reason the calculated time histories 
sometimes anticipate the time of shock front arrival (e.g. 

Fig. 24.4, position 4.0, 0.4) and do not resolve weaker shocks 
subsequent to the first, although these may be detected occasionally 
as a rounded bump in the normally exponentially decaying curve. 
Efforts are being made to improve the space and time 

resolution of density and pressure calculated from the 

particle trajectories. 

The lack of resolution close to the shock front does not 
occur in the case of particle velocity, which can be measured 
with reasonable accuracy as soon as the shock has traversed 
the relatively small space represented by an individual smoke 
puff. This improved resolution is manifested also in the 
dynamic pressure histories which depend on particle velocity 


squared. 
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3.5 Comparisons with gauge results 


The hydrostatic and total pressure time-histories were 
calculated at gauge positions within the smoke puff grid and 
the results from the particle trajectory analysis are compared 
with the corresponding transducer outputs in Figures 27.1 
to 27.4. 

The agreement between the results from the two measurement 
methods is good although, as previously discussed, the poor 
time resolution of the particle trajectory results does not 
permit Sudeqed Peeaeies of multiple shocks. This is illustrated 
in Figure 27.1 for the 60.20 location. However, the agreement 
between the two impulse curves is excellent. 

The good agreement between the total-pressure time- 
histories from the two measurement methods shown in Figures 27.3 
and 27.4 further confirms the validity of the technique used 
to calculate the total pressure, as described in the addendum 


to Volumes 1 and 2 of this report. 


In considering the above comparisons it must be remembered 
that determination of pressure time-histories was not an 
original objective of the particle trajectory analysis project, 
but the reasonable agreement with gauge measurements gives 
some indication of the reliability of the method. Also the 
gauge measurements were made on the opposite side of the 
charges to the smoke puff grid so that some differences might 


be expected due to slight non-symmetries of the blast waves. 
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Although the effect is more clearly seen from the gauge 
results rather than the particle trajectory analysis, it is 
interesting to note the difference in the blast wave signature 
at location 60.10, 10ft above the rough ground, compared with 
that at 60.40, 10ft below the interaction plane (Fig. 27.1). 
Above the ground the two initial shocks indicate that the 
gauge was above the triple point and detected the primary and 
the reflected shock. Below the interaction plane the single 
shock indicates that the gauge was in the Mach stem below 
the triple point. This confirms the photogrammetrical measure- 
ments (Dewey, et al., 1975) which showed the triple point above 
the ground at this location to be at height of approximately 
5.8 ft and to be approximately 9.6 ft below the interaction 


plane, i.e., almost coinciding with the gauge position. 


3.6 Comparisons of time histories at different heights 


Dewey et al., (1975) pointed out that although the 
shock strength below the interaction plane at a given distance 
from the vertical charge axis was greater than it was above 
the rough ground, the gauge results showed that the hydrostatic 
pressure impulse was greater along the ground. In other 
words, along the ground the peak pressure value at the shock 
was less, but the decay of the pressure was slower and the 


duration longer. 
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The results from the particle trajectory analysis have 
been studied to see if these conclusions are confirmed. 
Unfortunately, gaps in the smoke puff grid close to the 
roughened ground, seen in figure 22, meant that fewer results 
were obtained in this region than had been hoped for. Neverthe- 
less, a number of comparisons between time histories below the 
interaction plane and above the ground are possible. The 
appropriate time histories to be compared are the particle 
velocity profiles in figure 23.2 at (2.5, 1.6) and (2.5, 0.2), 
in figure 23.3 at (3.0, 1.6) and (3.0, 0.2), and in figure 
23.4 at (4.0, 1.6) and (4.0, 0.3); the density profiles in 
figure 24.3 at (3.0, 1.6) and (3.0, 0.4) and in figure 24.4 at 
(4.0, 1.6) and (4.0, 0.4); the hydrostatic pressure profiles 
in figure 25.3 at (3.0, 1.6) and (3.0, 0.4), and in figure 
25.4 at (4.0, 1.6) and (4.0, 0.4), and the dynamic pressure 
profiles in figure 26.3 at (3.0, 1.6) and (3.0, 0.4) and in 
figure 26.4 at (4.0, 1.6) and (4.0, 0.4). Although the 
differences are not great when the accuracy of these profiles 
is considered, it appears to be significant that in every one 
of the nine comparisons, close to the rough ground the shock 
front arrives later and the peak value is less than below the 
interaction plane, but the compared history profiles always 
cross so that in the latter part of the wave both the hydrostatic 
and the dynamic pressure are greater above the ground than at 


the interaction plane. 
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